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Abstract 

Non-linear recurrences which generate integers in a surprising way 
have been studied by many people. Typically people study recurrences 
that are linear in the highest order term. In this paper I consider 
what happens when the recurrence is not linear in the highest order 
term. In this case we no longer produce a unique sequence, but we 
sometimes have surprising results. If the highest order term is raised 
to the m*^ power we expect answers to have m*^ roots, but for some 
specific recurrences it happens that we generate rational numbers ad 
infinitum. I will give a general example in the case of a first order 
recurrence with m = 2, and a more specific example that is order 3 
with m = 2 which comes from a generalized Somos recurrence. 

1 Introduction 

Many people have studied non-linear recurrences that generate sequences 
of integers despite the fact that every iteration of the recurrence requires 
division by some previous term in the sequence. These types of non-linear 
recurrences generally have the following form 

a{n) = L{a{n — 1), ... , a{n — k)) (1) 

where L is a Laurent polynomial with integer coefficients, i.e.- L is in the set 
Z[xf^, . . . ,a;^^]. Well studied examples of this phenomenon are the Somos 
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sequences, introduced by Michael Somos in 1989 [2], defined by the recurrence 




i=l 



with initial conditions s(m) = 1 for m < k. For k = 2,3 the recurrence 
generates the infinite sequence {l}^^]^. More interestingly, for = 4, 5, 6, 7 it 
is known that these recurrences each generate an infinite sequence of integers 
([1][2])- There are, of course, other examples of this integrality phenomenon, 
and many are generalizations of the Somos recurrences (for some examples 
see [2]). 

It seems to be the case that all recurrences of the form that have been 
studied have no exponent on a{n). In this paper I will discuss recurrences of 
the form 



where m > 1, and L is still a Laurent polynomial. I will refer to m many times 
in this paper, and (unless otherwise stated) this will refer to the exponent 
on a{n) in the left-hand side of a recurrence of the form ([2]). 

In general one would imagine that if a{n) is raised to a power m > 1 then 
what is generated is not a sequence at all, since solving an equation of degree 
m yields up to m answers. So I need to introduce the concept of a recurrence 
tree. 

Definition 1. A recurrence tree is a way of storing the values generated by 
a recurrence of the form ^ with m > 1. Solving for the n + 1** term, given 
a specific n*^ term requires solving an equation of degree m. This yields up to 
m possibilities for the n + P* term. We can store these numbers in a complete 
m-ary tree. 

For example when m = 2 we get the following structure 



a(n) 



m 



L{a{n), a{n — 1), ... , a{n — k)) 



(2) 



a(l) 



a(2)i a(2)2 

/ \ / \ 

a(3)i,i a(3)i,2 a(3)2,i 0(8)2,2 
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Also, since solving an equation of degree m yields answers which can be 
in C, we may expect that the numbers generated are not rational. However, 
in some cases a recurrence of this form generates rational numbers. When 
the tree consists only of numbers in Q (resp. Z) we will call it a rational 
(resp. integer) recurrence tree. 

One way to come up with recurrences that obviously generate rational 
recurrence trees is to take a recurrence that generates integers and find the 
"ratios of ratios" sequence. 

Definition 2. Given the sequence {b{n)}'^^^, we call | ^^b(n)^ | 

quence of ratios of {b{n)} and | ''^b(n+i)/b(n)^ } sequence of ratios of 

ratios of {b{n)}. 

Obviously, if a sequence {6(77,)}^^ C Z then the sequence of ratios of {b{n)} 
and the sequence of ratios of ratios of {b{n)} are in Q. 

Of course, it may not be the case that the recurrence that generates these 
ratio sequences has the m > 1 property, but in the case of the generalized 
Somos-4 sequences we can find an alternate recurrence for the sequence of 
ratios of ratios that does have this property. We can then generalize and find 
new recurrences that do not obviously generate rational numbers. 



2 Generalized Somos-4 Ratios of Ratios Se- 
quence 

Let {sc(n)}J^^ be a sequence defined by the following recurrence: 

Sc{n)sc{n - 4) = CiSc(n - l)sc{n - 3) + C2Sc{n - 2)^ (3) 

with initial conditions Sc{i) = 1 for i < 4, where c = (ci, C2) G . This is a 
special case of the three term Gale- Robinson recurrence ([I] [2]) that further 
specializes to the Somos-4 recurrence when ci = C2 = 1. The first few terms 
of the sequence are 

1, 1, 1, 1, Ci + C2, cl + C1C2 + C2, cl + 2clc2 + C1C2 + 2CiC2 + C2, . . . 

Using cluster algebras and the Caterpillar Lemma, Fomin and Zelevinsky 
proved that the recurrence Q generates a sequence of integers pQ. 
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Now, define sequences {tcin)}'^^^, and {ac(^)}$^i by 

tc{n) = Sc{n + l)/sc(n) 
ac{n) = tc{n + l)/tc{n) 

then {tc(^)}^i is the sequence of ratios of Sc{n) 

cj + C1C2 + C2 + 2clc2 + C1C2 + 2cicl + cl 



{tcin)} = <^ 1, 1, l,ci + C2, , 2 , 

t Ci + C2 Cj + C1C2 + C2 

and {ac(n)} is the sequence of ratios of ratios of Sc(n). 

r / \Too fi 1 C? + C1C2 + C2 1 

{ac(^)}„=i = <l,l,ci + C2,—^——^,...> 

In this paper we will be interested in the sequence {ac{n)}. By algebraic 
manipulation we can easily find a first order quadratic recurrence for adn). 

Claim 1. The sequence {ac{n)}'^^i is defined by the recurrence 

ttcin + 2)ac(n + l)^ac(n) = Ciadn + 1) + C2 (4) 

with initial conditions ac(l) = ac(2) = 1. 

Proof. We will simply manipulate the recurrence equation for Sc{n) to look 
like the recurrence equation (jl]). 

Scin + A)sc{n) = CiSc{n + 3)sc{n + 1) + C2Sc{n + 2)^ 
Scjn + 4:)sc{n) _ Scjn + 3)sc{n + 1) 
s,(n + 2)2 """^ Se(n + 2)2 

Notice that the Sc term on the right side is adn + 1). By multiplying and 
dividing by the correct terms on the left side we will get the left side of @. 

Sc{n + A)sc{n) _ Sc{n + 4)sc(n) Sc{n + 2)'^Sc{n + 3)^Sc(n + 1)^ 
Sc(n + 2)2 ~ Sc{n + 2)2 Sc{n + 2ysc{n + 3)2sc(n + l)^ 

_ Scin + A)sc{n + 2) Sc(n + 3)Sc(^ + 1)^ + 2)sc(n) 
~ Sc(n + 3)2 Sc(n + 2)4 Sc(n + 1)2 

= ac(r;, + 2)ac(ra + l)^ac(n) 

Finally we obtain 

ac{n + 2)ac{n + l)^ac(n) = Ciac{n + 1) + C2 
which is (HI). □ 
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Unfortunately, this recurrence for the ratios of ratios of Sc{n) does not 
satisfy m > 1. The proof of the next claim will help to create a recurrence 
with m = 2. 

Claim 2. The sequence generated by the recurrence ^ also satisfies the 
recurrence 

ac{n + 2)ac{n + 1)^ + ac{n + l)^ac(n) = (2ci + C2 + l)ac{n + 1) — Ci (5) 

Proof. Showing the converse, that the sequence defined by ([5]) satisfies @, 
will prove this claim because of uniqueness of the sequence. So assume the 
sequence {ac(^)}5^Lo is defined by Now let us define a function T{n) by 

T{n) := ac{n+lYac{nY — {2ci+C2 + l)ac{n+l)ac{n) +Ciac{n+1) +Ciac{n) +02 

I claim that it is enough to show T{n) = 0, for if this is true then rearranging 
terms and dividing both sides by ac{n + l)^ac(n) we get 



Ciac{n + 1) + C2 = (2ci + C2 + l)ac(n + l)ac{n) — ciadn) — ac{n + lYadnY 
(2ci + C2 + l)ac{n + 1) — Ci — adn + l)^ac(n) 



ciac{n + 1) + C2 _ , ^ , iN„ ^„ , i\ „ „ ^„ , i\2 



Ciadn + 1) + C2 _ (2ci + C2 + l)ac(^ + 1) - Ci - adn + l)^ac(n) 
ac(r;, + lYadn) adn + 1)^ 

The RHS of the above equality equals aciji + 2) because we assumed the 
sequence {acin)} is defined by the recurrence ([5]). Therefore we also have 
that adn + 2) = LHS, i.e.- 

/ , ON ciadn + 1) + C2 
ac{n + 2) = — — 

ac{n + iyac{n) 

which is recurrence @. So the sequence {ac(n)}^^, generated by ([5]), also 
satisfies the recurrence (jl]). 

All that is left is showing, by induction, that T(n) = for all n. For 
n = 1 we do the following calculation 

r(l) = ae(2)V(l)' - (2ci + 02 + l)ae(2)ae(l) + cia,(2) + ciae(l) + C2 

= 1 ■ 1 - (2Ci + C2 + 1) • 1 ■ 1 + Ci ■ 1 + Ci ■ 1 + C2 
= 1 - (2Ci + C2 + 1) + 2Ci + C2 

= 
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Now assume that T{n — 1) = for some n. We substitute for adn + 1) in 
T{n) from and simplify to obtain 

T{n) =ac{n — l)'^ac{nY — (2ci + C2 + l)ac{n — l)ac(n) + Ciadn — 1) + 
+ Ciadn) + C2 
=T{n - 1) = 

Therefore by induction, T{n) = for all n and the claim is proved. This 
proof used ideas from Guoce Xin's paper [3]. □ 

Coming out of the proof of Claim [2] we get that T{n) = is a first order 
recurrence with m = 2 for {ac{n)} as we had hoped. One would expect 
that, since {adn)} is by definition a single sequence, the recurrence tree for 
T(n) = would somehow consist only of this single sequence. Indeed this is 
the case as we prove now. 

Claim 3. The recurrence tree for 



with a(l) = 1, produces a single sequence in the sense that at every level of 
the tree there is only one value that we haven't yet seen. 

Proof. Let X := ac{n + 1) and Y := ac{n), then the first order quadratic 
recurrence for the generalized Somos-4 sequence is rewritten as 



Given some value yo for Y there are two possible values for X which satisfy 
([8]). This corresponds to the fact that given some a^in) there are two possible 
values for ac{n + l). Using the quadratic formula, these two values, in terms 
of yo, are 



-(ci - (2ci + 02 + l)yo) ± v/(ci - (2ci + 02 + l)yo)' - Ay^iciyo + C2) 



Now, since these are values for a^n + 1) we substitute them back in for Y 
in (IH]), solve for X, and get potentially 4 possible values for ac{n + 2) that 



ac{n + l)^ac(^)^— (2ci + 02 + l)ac{n + l)ac{n) + 

+ Oia^n + 1) + Oiadn) + 02 = 



(6) 
(7) 



Y^X^ + (ci - (2ci + 02 + l)Y)X + (ciF + C2) = 



(8) 



6 



come from this specific adn) = yo- However, when we solve the quadratic 
equation 

(y^yx^ + (ci - (2ci + C2 + l)y+)X + (ci|/+ + C2) = 

for X the two solutions we get are yo and a large expression in terms of 
yo,ci,C2 (similarly for y~). This means that in the i^^ level of the tree, 
representing all possible values for ac{i), from each of the terms in the i — 1^* 
level there is at most one term that we haven't yet seen. Now, lets look at 
the second level given the initial condition (the root) ac(l) = 1. We solve 
the quadratic equation 

l^X^ + (ci - (2ci + C2 + 1) • 1)X + (ci ■ 1 + C2) = 
X^ + (-ci - C2 - l)X + (ci + C2) = 
X = 1 or ci + C2 

So on the second level we only have one new term. Therefore, on the third 
level, and all subsequent levels, we also only have one new term. □ 

Since the recurrence tree for ([6]) consists only of numbers from the se- 
quence of ratios of ratios of {sc{n)}, it must be a rational recurrence tree. 
So we have found an example of a non-linear recurrence with m > 1 that 
generates rational recurrence tree. However, this example was constructed 
in such a way that it had to generate a rational recurrence tree. In the next 
section I will generalize this example to get nontrivial sequences generating 
rational recurrence trees. 



3 Generalized First Order Quadratic Recur- 
rence Tree 

The general form of a first order non-linear recurrence is 

m 

Y^Pi{a{n))a{n + iy = (9) 

i=0 

where Pi{Y) is a polynomial in Y of some degree di. For example, the 
sequence of ratios of ratios of generalized Somos-4 has recurrence given by 
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(Q where m = 2 and 



P2{Y) = 

Pi(F) = ci-(2ci + C2 + l)F 
Po{Y) = C2 + c,Y 



For the remainder of this section we will assume that m = 2, do = di = 1, 



where Ai, A2, Bi, B2 G Z. Under certain minimal sufficient conditions a 
recurrence of this form will generate a rational recurrence tree. 

Proposition 4. Let a(l) = 1 in the recurrence ^ with coefficient polyno- 
mials given by [T^) . If 

(i) Ai = B2 and 

(a) solving for a{2) yields rational numbers, 

then the recurrence generates a rational recurrence tree. 

Proof. First we will show that Ai = B2 implies that for every term ai = a{n) 
coming from solving 

a{nfa{n - if + {Ai + A2a{n))a{n - 1) + {Bi + Aia{n)) = 

with a{n — 1) = Oq, we get only one new a{n + 1). In other words, solving 

a{n + ifaj + {Ai + A2a{n + l))ai + {Bi + Aia{n + 1)) = 

for a{n + 1) yields the solutions {a{n + 1), ag}. In recurrence tree form this 
looks like: 



and P2{Y) = Y\ Let 



P2{Y) = Y^ 
PiiY) = Ai + A2Y 
PoiY) = B, + B2Y 



(10) 



a(n — 1) = oq 



/ 

a{n)i = ai 

/ ^ 




a(n + 1)1,2 = ao + l)i,i 
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Let Oo be a term in the n — 1** level. Its children in the recurrence tree 
are the solutions of the quadratic equation 



alX^ + {Ai + A2ao)X + {Bi + AiUo) = 



In other words, possibilities for a{n) given that a{n — 1) — clq are 



a(n)i 



2al 



a{n)2 



{Ai + A2Q0) - v/(Ai + AaQp)^ - 4ag(Ei + Aia») 

2a2 



If a{n) = a{n)i, for i = 1,2, then to get a(n + 1) we solve 

a{n)^a{n + if + {Ai + A2a{n),)a{n + 1) + {Bi + Aia{n)i) = 

for a{n + 1). The goal is to show that is a solution for a{n + 1), so it is 
enough to show that a{n)^al + {Ai + A2a{n)i)ao + (-Bi + Aia{n)i) = for 
i = 1,2. This is nothing but algebraic manipulation that can be easily done 
using Maple or any other computer algebra system. 

The other assumption that we made is that both answers for a(2) are 
rational. When we solve for level 3 we know that in each case one answer 
will be from level 1, so rational. Then the other answer must be rational 
because the product of the two answers is the constant term in the quadratic 
polynomial, Bi + y4ia(2)j, which is rational. Likewise, if levels n — 1 and n 
are rational then level n + 1 will be rational. So by induction we see that all 
levels are rational numbers. □ 

Even with the stipulation that Ai = B2 and a(2)i, 0(2)2 are both ratio- 
nal, this more general first order recurrence encompasses more than just the 
Somos-4 ratios of ratios. For example, let Ai = i?2 = 1, ^2 = 5, -Bi = 8 with 
the initial condition a(l) = 1, then the recurrence is 



a(n - lfa{nf + (1 + 5a{n - l))a{n) + (8 + a{n - 1)) = (11) 



The reader may check that a(2)i, 0(2)2 are in fact rational. There is no 
(ci,C2) such that IQ, the generalized Somos-4 ratios of ratios recurrence, is 
the same as ffTTj) . 

Another fact that is worth pointing out is that these recurrences produce 
at most two sequences. In Section [2] we showed that the first order quadratic 



9 



recurrence for ac{n), the ratios of ratios of the generahzed Somos-4 sequence, 
generates a unique sequence. This was partly due to the fact that at level 2 
one of the solutions is ac(l), so we only have one new solution. However, it is 
possible to have two new solutions on level 2. In this case, on level n > 2 we 
could have two new solutions, each coming from one of the a(2). However, 
we can never get more than two sequences. 

4 Recurrence Tree for Generalized Somos-4 
Sequence 

So far we have looked at first order nonlinear recurrences that generate trees. 
What about higher order nonlinear recurrences that generate trees? As an 
example we will look at the generalized Somos-4 recurrence (j3]). This re- 
currence is order 4 and quadratic, but with m = 1. Since we have a first 
order quadratic recurrence for the ratios of ratios of the generalized Somos-4 
sequence given by we can "unfold" this recurrence to get an order 3 re- 
currence that should be satisfied by the generalized Somos-4 sequence. Since 
we assumed that ac(n) := tc{n+l)/tc{n) and tc{n) := Sc(n+ l)/sc(n) we can 
substitute into (Q. The recurrence we obtain is 



-((2ci + C2 + l)sc(n + 2)sc(n + l)sc(n) - CiSc{n + l)^)sc(n + 3) = 

In Section [2] we were able to show that the recurrence tree for the ratios of 
ratios of the Somos-4 recurrence is just a sequence in disguise because there 
is only one new term per level. One might think that the same happens here 
since f lT^ is supposed to generate the generalized Somos-4 sequence, however 
this is not the case. Instead, we get many new terms per level, and therefore 
generate many sequences. Also surprisingly, we even get some non-integer 
rational numbers in the tree. As it turns out, at least one of the sequences 
in the tree has a simple closed form. 

Proposition 5. One of the sequences generated by the recurrence fT^) is 



Sc{n + 3fscinf + CiSc{n + 2)^Scin) + C2Scin + 2)^Sc(ra + 1)^ 



(12) 



s,(l) = s,{2) = s,(3) 



l,{(ci + C2)^(")} 



n=4 



oo 



where f{n) 



(n-3)^ 



4 
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Proof. First note that the sequence {f{n)}'^^^ is 

0,1,2,4,6,9,12,16,20 

I claim that f{2i + 1) - /(2i) = i - I, and /(2(i + 1)) - f{2i + 1) 
The vakie of f{2i) is 



i - 1. 



(2z - 3)^ 



and the value of f{2i + 1) is 

:2z + 1 - 3)2 



r - 3z + 2 



2i + 1 



Their difference is i — 1 as claimed. Now, /(2(i + 1)) is 

4 J 

The difference f{2{i + 1)) — f{2i + 1) is indeed i — 1. 

Now we will prove this proposition by induction. First, the base case: we 
need to show that Sc(4) = (ci + C2)°,Sc(5) = (ci + C2Y,Sc{6) = (ci + 02)^. 
This is straightforward computation, substitute in Sc(l) = Sc{2) = Sc(3) = 1 
in (fT2i) and solve for Sc(4). We get two possibilities, 1 and ci + C2. In this 
case we choose Sc(4) = 1. Now we again substitute Sc{2) = Sc(3) = Sc(4) = 1 
in f fT2l) and solve for Sc(5). This is the same quadratic equation so we get 
the same 2 solutions, this time we choose Sc(5) = Ci + C2. To finish off the 
base case we substitute Sc(3) = Sc(4) = 1, ■Sc(5) = Ci + C2 in f[T^ and solve 
for Sc(6). Our two possible solutions are 

cj + C2C1 + C2 and (ci + 02)^ 

so we choose Sc(6) = (ci +02)^. So the base case is true. Now assume, as the 
inductive hypothesis, that the proposition is true up to n + 2. We have two 
possibilities 



Case 1: n = 2i 



Sc{2l) =(C1 + C2f 

S,{2l + l)={c^+C2f+'-^ 
S,(2(2 + l))=(Ci+C2)^+2.-2 
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If we substitute this into (fT^ and simplify, we obtain the quadratic 
equation 

s{2t + 3)2(ci + - ((ci + C2 + l)(ci + C2f'~'''-')s{2t + 3) + 

+ (ci + C2)(c2 + Ci)^'=+6^-^ = 

which we can easily solve using the quadratic formula. Our two possi- 
bilities when solving for s{2i + 3) are 

+ 3) = I + l'^k+3^-3 

We expected Sc(2(i + 1) + 1) = (ci + 02)''+^'-'^+' = (ci + 02)^'+^*"^ and 
we have it if we choose the "+" in the quadratic formula. 

ase 2: n = 2i + 1 

S,{2i + 1) =(ci + C2)*^ 
Se(2(z + l))=(ci + C2)''+^-^ 
Se(2(z + l) + l)=(ci + C2)'=+'*-^ 

Again we substitute this in to f[T^ to obtain the following quadratic 
equation 

s{2i + 4)2(ci + C2)''+ 

- ((2ci + C2 + l)(ci + c2)3''+=^*-' - ci(ci + 02)=^'=+=^'"') s(2z + 4) + 

+ Ci(Ci + 02)^'+^*-=^ + C2(C2 + Ci)^'=+^*-^ = 

which we can again solve using the quadratic equation and obtain 

S[ZZ + 4j - I ^ C2)'=+3^-3(c2 + C2C1 + 4) 

Again, if we choose the " + " in the quadratic equation we get 

s,{2{z + 2)) = (ci + c2)'=+2^-i+* = (Ci + C2)'+^'-' 
as expected. 

□ 
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The fact that we were able to find a nice closed form for one of the integer 
sequences in this recurrence tree is very surprising. The closed form for the 
generalized Somos-4 sequence is in terms of elliptic theta functions [3], but 
by finding a lower order recurrence with higher degree we were able to find 
a polynomial sequence. 

In looking at the tree that the recurrence f[T^ generates I have noticed 
that something even more general seems to be true. 

Conjecture 6. Let T be the recurrence tree for [T^) . Every integer, and the 
numerator and denominator of every (reduced) non-integer rational number 
in T , are products of terms in the generalized Somos-4 sequence. 

Clearly, proposition [S] is consistent with this conjecture since Sc(5) = Ci + C2 
in the general Somos-4 recurrence ([3]). 

5 Maple code 

This subject could not have been studied without the use of a computer alge- 
bra system, Maple in my case. The Maple code accompanying this paper can 
be found on my website http://math.rutgers.edu/~eahogan/maple/. I 
created programs that calculate the recurrence tree for a given recurrence 
of any order and any degree. These programs can be found in the file 
RecurrenceTree.txt. I also created programs that generate all recurrences 
that have rational recurrence trees to a certain depth (i.e.- if the tree for a spe- 
cific recurrence is rational up to a test depth, the program outputs that recur- 
rence). Those can be found in file GenerateRationalRecurrenceTrees . txt. 

6 Conclusion 

My study of recurrence trees in this paper may be just the beginning. The 
first order recurrences I looked at were limited to the case where m = 2. 
As m grows it seems less likely that the recurrence trees generated will be 
rational. Though it could be the case that a subtree is rational, or perhaps 
just a single sequence. The only higher order recurrence tree I investigated 
was that of the generalized Somos-4 recurrence. That specific recurrence is 
not completely characterized, but I suspect a generalization, along the lines 
of section [3l can be made which may yield behavior like (fT2!) . 
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